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ABSTRACT 

We have discovered small whirlpools in the Sun, with a size similar to the terrestrial hurricanes 
(< 0.5 Mm). The theory of solar convection predicts them, but they had remained elusive so far. The 
vortex flows are created at the downdrafts where the plasma returns to the solar interior after cooling 
down, and we detect them because some magnetic bright points (BPs) follow a logarithmic spiral in 
their way to be engulfed by a downdraft. Our disk center observations show 0.9x10^^ vortexes per 
Mm^, with a lifetime of the order of 5 min, and with no preferred sense of rotation. They are not evenly 
spread out over the surface, but they seem to trace the supergranulation and the mcsogranulation. 
These observed properties are strongly biased by our type of measurement, unable to detect vortexes 
except when they are engulfing magnetic BPs. 

Subject headings: convection - Sun: photosphere - Sun: granulation 



1. INTRODUCTION 

A proper understanding of the solar surface convection 
has impact on many different areas of astrophysics. For 
example, (a) the solar metallicity is used as a univer- 
sal reference, and it has been recently modified by fifty 
percent after considering re alistic convective motions in 
the analysis (|Asplundll2005l ). and (b) turbulent dynamo 
action can produce ma gnetic fields in many astrophys- 
ical e nvironments (e.g.. iMaron et al]l2004t (Bercik et al.l 
|2005( ). but its observational study is so far confined to the 
turb ulent dynamo dr i ven by the solar convectiv e motions 
re.g.. lCattaneol [19991 : iVogler fc Schiissleil[200l . 

We have discovered small whirlpools in the Sun, with 
a size similar to the terrestrial hurricanes. This dis- 
covery confirms a specific prediction of the current the- 
ory of solar convection that had remained elusive so 
far. The energy produced in the solar interior is trans- 
ported to the surface firstly by radiation exchange and, 
during the last t hirty percent of the way, by convec- 
tive motions (e.g.. lSti:>3 ll991). The hot buoyant plasma 
rises, releases energy, and then falls down. Accord- 
ing to the current theory, these convective motions are 
driven by st rong highly l ocalized downdrafts at the so- 
lar surface ('S pruit et all Il99(]t IStein fc Nordlundl 1 19981 : 
lRasti ri998). The downdrafts are sinks where the cold 
plasma returns to the solar interior. Since the mat- 
ter has angular momentum with respect to the drain- 
ing point, it must spin up when approaching the sink, 
giving rise to a whirl flow (bathtub effect). Although 
the vortex motions predicted by the theory have been 
repeatedly soug,ht, they had not been found so far 
(lNovembeijll98 9': 'Spruit et al.lll99(]l: iRoudier et al.l ll997t 
iHoekzema eTal . 1998; Nis enson et al.ll2003D . Tflie ex- 
ception in iBrandt et al.l 1198^ represents the same phe- 
nomenon at much larger scale.) We have detected 
these vortexes as proper motions of magnetic bright 



points (BPs), which follow spiral paths in the way 
to be engulfed b y a downdraft. Th e BPs are easily 
blurred by seeing (|Title fc BergeHfl996l) . thus the unique 
90 km spatial resolution provided by the Swedish So- 
lar Telescope at La P alma turned out to be necessary 
(jScharmer et al.ir2002t l. Magnetic BPs are present almost 
every\ Yhere on the Sun, also in the quiet internetwork re- 
gions ([Sanchez Almeida et al.]|2004l : Ide Wiin et al.ll2005l 

l200 8f). Precisely looking for properties of these quiet 
Sun BPs, we found out that some of them spiral when 
heading towards the supposed location of a sink, which 
is the result reported in this Letter. We note that 
the vortical motion of photospheric magnetic concentra- 
tions may have significant impa ct on its own, e. g., on 
the heating of the solar corona ([van Ballcgooiicn et al.l 

119981: iGudiksen fc Nordlundl 120051 : lAmari et al.l i2008). 
Previous searches for phot ospheric vortexes were often 
motivated by this iii terest (jvan Ballegooijen et al.lll998l : 

iNisenson eralll2003[ ). 

2. OBSERVATIONS 

A quiet Sun internetwork target was observed near so- 
lar disk center with t he Swedish Solar Telescope (SST; 
IScharmer et al.l l2002[ ) . During observation seeing was 
fluctuating with minutes of excellence (September 29, 
2007; from 9:13:57 to 9:45:55 UT). We used 12 bit Ko- 
dak Mega Plus II 2048 x 2048 cameras, which render an 
image scale of 0.034 " pix^i. The field-of-view (FOV) 
includes various supergranules. We observe simultane- 
ously in the G-band (a 10.8 A wide filter centered at 
4305.6 A), in Ca H line-core, and in the continuum close 
to Ca H. The images were reconstructed with multi-frame 
bhnd deconvolution (|van Noort et al.ll2005D to produce 
movies whose snapshots have an angular resolution close 
to the diffraction limit of the SST at the working wave- 
length (^O'.'l). Each snapshot results from combining 
J.25 images, rendering a mean cadence of one snapshot 
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Fig. 1. — Sequence of snapshots corresponding to the event ana- 
lyzed in detail in the Letter. The snapshots are separated by 15 s, 
and the bar on the lower left image corresponds to 1 Mm on the 
Sun. We have marked with symbols two sets of two nearby BPs, so 
that it is easy to appreciate how the two BPs rotate with respect 
to each other. (Although the sequence is continuous, the two up- 
per rows show only one pair, whereas the two lower rows show the 
other pair.) Note how the bright structures in each snapshot depict 
a spiral pattern resembling the vortex derived from BP proper mo- 
tions. The animation given as on-line material allows to appreciate 
the vortex motions in more detail. 

3. WHIRLPOOLS 

The G-band movie often shows BPs swirling around 
intergranular points where several dark lanes converge. 
These motions are reminiscent of the bathtub vortex 
flows predicted by numerical simulations of convection, 
and which are driven by the granulation downdrafts (§[!]). 
Since we have no other information but the proper mo- 
tions of the BPs, we cannot proceed without assuming 
that these proper motions trace true plasma motions. 
We cannot think of any other physical mechanism able 
to generate these very special trajectories, but it must be 
clear that treating the proper motions as plasma motions 
is a working hypothesis'^. 

We will discuss statistical properties of the whirlpools 
in § m This section is focused on a particularly well de- 
fined event whose snapshots are shown in FigHJ (See also 
the animation provided as on-line material.) Up to eight 
different BPs move under the influence of this whirlpool 
during the 5 min duration of the event. Their trajectories 
are shown in Fig. [51 where we use a different symbol for 
each trajectory, i.e., the set of positions of a single BP 

* One may argue that BPs are not good passive tracers because 
they posses strong kG magnetic fields and, therefore, the mag- 
netic forces back-react on the plasma flows counteracting the drag 
forces. This theoretical prejudice has been proved to be erroneous, 
since the BPs move exactly as dic tated by the granular flow field 
(see Ivan Ballegooiien et al.lll998l) . The reason seems to be the 
small cross section of the magnetic concentrations. Drag forces 
overcome magnetic forces when the magnetic concentration s are 
slender enough (e.g.. ISchiissler|[l986l : [Sdnchez Almeidall2001h . 



as time goes on. The positions have been determined 
as the centroid of the BP in successive snapshots. The 
identification of each particular BP has been carried out 
by eye, playing back and forth the movie with the FOV 
of the event. It is clear from Fig. [5^ that the BPs swirl 
around a quite stable central point. Only two BPs follow 
straight trajectories (see Fig.l^Jj), but this path coincides 
with the position of one of the intergranular lanes that 
feeds the whirlpool, and which channels the BP motion 
(see also the top row of Fig. [T]). 

If the analogy with a bathtub vortex holds for the ob- 
served whirlpools, then a passive tracer should approx- 
imately follow a logarithmic spiral when heading to the 
sinkhole, i.e., 

9-9o=AHr/ro), (1) 

where r and 9 stand for the radius and azimuth in polar 
coordinates with the draining point as origin. The sym- 
bol 9q corresponds to the azimuth at the arbitrary radius 
ro, and 

A = Ue/Ur, (2) 

is a constant of motion given by the ratio between the 
azimuthal velocity Ue and the radial velocity Ur- The pa- 
rameter A controls the shape of the spiral, which becomes 
a straight line for A equals zero, and a circle for A —^ oo. 
Equation ([T]) follows from the conservation of mass and 
angular momentum in a non- visco us incompressible fluid 
moving in a bathtub vortex (e.g.. iFevnman et al.lll96S 
lMassevlll983[ ). These conditions are only loosely met by 
the solar plasma, however, the observed motions match 
quite well this approximate behavior. Using a non-linear 
least squares procedure, we have fitted the trajectories 
of the BPs in Fig. [2^. We assume the whirlpool to be 
the same for all of them, i.e., all trajectories share a com- 
mon origin and the same A. Each BP trajectory has its 
own starting point parameterized by ^ at tq, i.e., by ^o- 
The fitting is fairly robust and independent of the initial 
conditions, providing, 

A = 6.4 ±1.6, (3) 

where the mean value corresponds to the fit in Fig. 
(the solid lines), and the error bar indicates the stan- 
dard deviation among the inferred values when each one 
of the BP trajectories is excluded from the fit. The area 
of the swirling region turns out to be of the order of 
0.5x0.5 Mm^. The size of the sinkhole (i.e., of the hypo- 
thetical downdraft) has to be smaller than the radius of 
curvature of the smallest trajectory, which corresponds 
to the one marked by triangles in Fig. (5] Its radius of 
curvature is of the order of 100 km, i.e., the downdraft 
seems to be spatially unresolved to SST observations. BP 
proper motion velocities span from one to four km 
and, with some noise, they tend to increase towards the 
sinkhole. 

4. STATISTICAL PROPERTIES 

Detained visual inspection of the G-band movie shows 
the whirlpool motion described above to be quite com- 
mon. In order to have a statistically significant descrip- 
tion of the phenomenon, we carried out a careful visual 
analysis of the full time series looking for swirl motions. 
Roughly speaking, one finds two types of motions that 
we have identified as vortexes: (1) two or more BPs move 
(in the same or in opposite directions) towards a point, 
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Fig. 2. — (a) Logarithmic spiral (the solid lines) that fits the 
trajectories of six observed BPs (the symbols), (b) Same as (a) 
but with a FOV twice larger and including two additional BPs (the 
solid circles) whose paths are not set by the whirlpool, but they 
follow a nearby intergranular lane. The square in (b) represents 
the FOV in (a). (The assignment of symbols to BPs differs from 
that used in Fig.[l]) 

where they rotate with respect to one another, and (2) 
the rotation also includes the diffuse bright borders of the 
sometimes fragmented granules surrounding the turning 
point. Taking this vague definition in mind, we find 138 
swirls (Fig. [3]), 83% of which belong to type 1, and the 
rest to type 2. (The motion shown in Fig. [1] belongs to 
the second type.) Our FOV covers 49.7Mmx49.7Mm. 
The time series lasts 31.5 min. This corresponds to an 
space-and-time density equals to c? ~ 1.8 x 10~^ vortex 
Mm~^ min"-^. 

The lifetimes of the vortexes were also estimated. The 
starting time was defined as the time where the BPs be- 
gin their swirling motions. These motions often end up 
producing a central still BP. Then we take the end time 
as the moment when this BP starts fading away or frag- 
menting. If the central BP is not formed, our final time 
is the time when circular motions are no longer distin- 
guished. The mean lifetime of the 138 whirlpools turns 
out tober~(5.1±2.1) min, with the error bar repre- 
senting the standard deviation among all single measure- 
ments. Assuming that the number of swirls is stationary, 
at any given time one finds rd whirlpools on the solar 
surface, which corresponds to 0.90x10"'^ vortex Mm~-^. 
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Fig. 3. — Positions of the observed vortexes on the average G- 
band image. They trace a pattern resembling the supergranulation 
and the mesogranulation. Plus signs correspond to counterclock- 
wise motions whereas circles correspond to clockwise motions. We 
do not think that such a distribution reveals an intrinsic property 
of the vortexes, but it is a bias set by using magnetic BPs to trace 
small scale proper motions. 



The lifetime estimates are only lower limits to the true 
lifetimes, since the vortexes should exist and last even if 
there are no magnetic BPs that allow us to detect them. 
For the same reason the number of detected vortexes also 
represents a lower limit. 

There is the same probability of clockwise and coun- 
terclockwise motions. We find 50.4% of the whirls with 
counterclockwise rotation, and 49.3% with clockwise ro- 
tation (Fig. [3]). 

The observed swirls are not uniformly distributed 
throughout the FOV. They follow a pattern that closely 
resembles supergranulation and mesogranulation cells 
(Fig. [3]). It is not clear whether this distribution re- 
veals an intrinsic property of the vortexes. Rather, 
it seems to be a selection effect we impose when us- 
ing G-band BPs as passive tracers of plasma mo- 
t ions. It is well known that BPs trace s upergranulation 
(Dunn & Zirker"1973"; 'Mehltrctter"'1974) and mesogran- 
ulation (Sanchez Almeida et al. 2004). Vortexes outside 
these areas lack of enough BPs and they will not be de- 
tected. 

5. DISCUSSION 

Observing these convectively driven vortex flows has 
been elusive because of the size of the whirlpools. Mo- 
tions take place at sub-arcsec scales and, therefore, close 
to the best angular resolution achieved at present. In 
addition, the case analyzed in §[3] seems to be extreme in 
the sense that the ratio between the azimuthal velocity 
and the radial velocity is particularly large {A ~ 6.4), 
which renders closed spiral paths and clear swirling mo- 
tions. Assuming that the velocity of the plasma that 
feeds the whirlpools has an important stochastic com- 
ponent set by the complex granular dynamics, then one 
should expect the azimuthal and the radial velocities to 
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Our work should be complemented with a spectro- 
scopic confirmation of the vortex motions. It requires 
observing the whirlpools out of the disk center so that 
the horizontal circular motions have a component along 
the line-of-sight (LOS). Such spectroscopic observation is 
challenging. Achieving the required spatial resolution is 
even more difficult in spectroscopy. Moreover, one would 
like to maximize the velocity component along the LOS 
by observing well outside the solar disk center. How- 
ever, the solar surface is corrugated and, outside the 
disk center, the granules protrude hiding the intergran- 
ular l anes that harbor t he vortexes (e.g.. ICarlsson et all 
[200l[KeIie"r et al.l[2?)0l . 



Fig. 4. — BP trajectories like the one analyzed in §[3] (the solid 
line, corresponding to A = 6.4), and what we expect to be the most 
probable case, with A ^ 1 (the dashed line). The swirling motions 
of the latter are difficult to detect since the deviations from straight 
lines are of the order of the angular resolution of the observation 
(represented by the hashed disk at the core of the whirlpool) . 

be similar, rendering ^ 1. In this case the swirling 
motions are moderate and so easy to overlook. The BP 
trajectories corresponding to these more typical cases are 
hard to distinguish from straight lines within the angular 
resolution of the observations (see Fig. |4]). 
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